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Abstract 

In this paper, we investigate a system of two nonlinear partial differential equations, arising from a 
model of cellular proliferation which describes the production of blood cells in the bone marrow. Due 
to cellular replication, the two partial differential equations exhibit a retardation of the maturation 
variable and a temporal delay depending on this maturity. We show that this model has a unique 
solution which is global under a classical Lipschitz condition. We also obtain the positivity of the 
solutions and the local and global stability of the trivial equilibrium. 

Keywords: nonlinear partial differential equation, age-maturity structured model, blood production sys- 
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1 Introduction 

We analyse, in this paper, a mathematical model arising from the blood production system. It is based 
on a system proposed by Mackey and Rudnicki [19], in 1994, to describe the dynamics of hematopoietic 
stem cells in the bone marrow. The origin of this system is a model of Burns and Tannock [2| (1970) in 
which each cell can be either in a proliferating phase or in a resting phase (also called Go-phase). The 
resulting model is a time-age-maturity structured system. 

Proliferating cells are in the cell cycle, in that they are committed to divide at the end of the mitosis, the 
so-called point of cytokinesis. After division, they give birth to two daughter cells which enter immediatly 
the resting phase. Proliferating cells can also die by apoptosis, a programmed cell death. 

The resting phase is a quiescent stage in the cellular development. Cells in this phase can not divide: 
they mature and, provided they do not die, they enter the proliferating phase and complete the cycle. 

The model in [H] has been analysed by Mackey and Rey [HIIIH] in 1995, Crabb et al. 011] in 1996, 
Dyson et al. [10] in 1996 and Adimy and Pujo-Menjouet [S] H] in 2001 and 2003. In these studies, the 
authors assumed that all cells divide exactly at the same age. 
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However, in the most general situation in a cellular population, it is believed that the time required for 
a cell to divide is not identical between cells (see Bradford et al. [B]). For example, pluripotent stem cells 
(which are the less mature cells) divide faster than committed stem cells, which are the more mature stem 
cells. In 1993, Mackey and Rey [H] considered a model in which the time required for a cell to divide 
is distributed according to a density, but the authors only made a numerical analysis of their model. 
Dyson et al. [HHH], in 2000, considered a time-age-maturity structured equation in which all cells do 
not divide at the same age. They presented the basic theory of existence, uniqueness and properties of 
the solution operator. However, in their model, they considered only one phase (the proliferating one), 
and the intermediary flux between the two phases is not represented. In 2003, Adimy and Crauste [2] 
considered a model in which the proliferating phase duration is distributed according to a density with 
compact support. They obtained global stability results for their model. 

In this work, we consider the situation when the age at cytokinesis depends on the maturity of the cell 
at the point of commitment, that means when it enters the proliferating phase. We assume that each cell 
entering the proliferating phase with a maturity m divides at age r = t(to), depending on this maturity. 
This hypothesis can be found, for example, in Mitchison [5T] (1971) and John [T^ (1981). This yields 
to the boundary condition (fTTj) . To our knowledge, nobody has studied this model, except Adimy and 
Pujo-Menjouet in [5], where they considered only a linear case. 

We obtain a system of first order partial differential equations with a time delay depending on the 
maturity and a retardation of the maturation variable. We investigate the basic theory of existence, 
uniqueness, positivity and stability of the solutions of our model. 

The paper is organised as follows. In Section [21 we present the time- age-maturity structured model. 
By using the characteristics method, we reduce this model to a time-maturity structured system, which is 
formed by two partial differential equations with a time delay depending on the maturity and a nonlocal 
dependence in the maturity variable. In Section [31 we first give an integrated formulation of our model by 
using the classical variation of constant formula and then we prove local existence of solutions, by using 
a fixed-point theorem, and their global continuation. We deduce the global existence. In Section IH we 
obtain the positivity of these solutions by developping a method described by Webb [24] . In Section [5l 
we concentrate on the stability of the trivial equilibrium of the system and, in the last section, we discuss 
the model and the asymptotic behaviour. 

2 Biological background and equations of the model 

Each cell is caracterised, in the two phases, by its age and its maturity. The maturity describes the 
development of the cell. It is the concentration of what composes a cell, such as proteins or other elements 
one can measure experimentally. The maturity is supposed to be a continuous variable and to range from 
m = to TO = 1 in the two phases. 

Cells enter the proliferating phase with age a ~ Q and they are committed to undergo cell division a 
time T later, so the age variable ranges from a = to a = t in the proliferating phase. We suppose that 
proliferating cells can be lost by apoptosis with a rate 7. 

At the cytokinesis age, a cell divides and gives two daughter cells, which enter immediatly the resting 
phase, with age a = 0. A cell can stay its entire life in the resting phase, so the age variable ranges from 
a = to a = -|-oo. The resting phase is a quiescent stage in the cellular development. In this phase, 
cells can either return to the proliferating phase at a rate f3 and complete the cycle or die at a rate 5 
before ending the cycle. According to a work of Sachs [55], we suppose that the maturation of a cell and 
the density of resting cells at a given maturity level determine the capacity of this cell for entering the 
proliferating phase. 

We denote by p(i, to, a) and n(t, m, a) respectively the population densities in the proliferating and the 
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resting phases at time i, with age a and maturity ni. The conservation equations are 

dp dp d{V{m)p) 



dt da dm 
dn dn d{V{m)n) 



--/im)p, 



= -(^S{m)+(3{m,N{t,m))'^n, (2) 



dt da dm 

where V{m) is the maturation velocity and N(t, m) is the density of resting cells at time t with a maturity 
level m, defined by 

p + OO 

N{t,m) — / n{t,m,a)da. 
Jo 

We suppose that the function V is continuously differentiable on [0,1], positive on (0,1] and satisfies 
V{0) = and 

ds 

/ TTTT = for ^ 6 (0' 1] • (3) 

Jo 

Since v^) ' '^ith mi < TO2, is the time required for a cell with maturity mi to reach the maturity m2, 
then Condition ([3]) means that a cell with very small maturity needs a long time to become mature. 
For example, if 

V{m) ~ am^, with a > and p > 1, 

then Condition ^ is satisfied. 

We suppose, throughout this paper, that 7 and S are continuous and non-negative on [0, 1]. The function 
(3 is supposed to be positive and continuous. 

Equations ^ and ^ are completed by boundary conditions which represent the cellular flux between 
the two phases. The first condition, 

r + co 

p(i,m,0)=/ P{m,N{t,m))n{t,m,a)da ^ P{m,N{t,m))N{t,m), (4) 
Jo 

describes the effiux of cells leaving the resting phase to the proliferating one. Cells entering the proliferating 
phase with age depend only on the population of the resting phase with a given maturity level. 

The second boundary condition determines the transfer of cells from the point of cytokinesis to the 
resting compartment. 

We assume that a cell entering the proliferating phase with a maturity m G [0, 1] divides at age T(m) > 0, 
and we require that r is a continuously differentiable and positive function on [0, 1] such that 

^'(™) + T77-T > 0' for e (0> !]■ (5) 
V (m) 

Since V{0) = 0, this condition is always satisfied in a neighborhood of the origin. If we suppose that the 
less mature cells divide faster than more mature cells, that is, if we assume, for example, that r is an 
increasing function, then Condition ^ is also satisfied. 

If one consider a cell in the proliferating phase at time t, with maturity m G (0,1], age a and initial 
maturity (that means at age a — 0) mo, then, naturally, we have 

mo < m and o. = — tt — ''"('7^o)• 

Jmo V{s) 

If m is the maturity of the cell at the cytokinesis point, then there exists a unique 0(m) G (0,m) (the 
maturity at the point of commitment) such that 

^-r(e(m)), (6) 
e(m) y[s) 
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because Condition ([5]) implies that the function 

ds 



V{s) 



is continuous and strictly decreasing from (0,r7i] into [— r(?7i), +00). Then, we can define a function 
e : (0, 1] (0, 1], where e(m) satisfies ©. 

From a biological point of view, Q(m) represents the initial maturity of proliferating cells that divide at 
maturity m (at the point of cytokinesis). Then, from the definition, the age of a cell with maturity m at 
the point of cytokinesis is t(9(to)). 

Remark that Q is continuously differentiable on (0, 1] and satisfies 

< 6(m) < m, for m E (0, 1]. 

This implies, in particular, that 

/■™ ds 

lim e(m) = and lim / — — = r(0) < +00. (7) 

rn^O ™^"Je(m) V(s) 

The property ([7]) means that cells with null maturity at the point of commitment keep a null maturity in 
the proliferating phase. 

The total number of proliferating cells at time t, with maturity m, is given by 

rr{e{m)) 

P{t,m) = / p{t,m,a)da. 



We consider the characteristic curves x : (— cx),0] x [0,1] [0,1], solutions of the ordinary differential 
equation 

{dx 
— {s,m) — V{x{s,m)), s < and m £ [0, 1], 
x(0,m) = m. 

They represent the evolution of cells maturity to reach a maturity m at time from a time s < 0. They 
satisfy x(s, 0) = and x{s, m) G (0, 1] for s < and m G (0, 1]. 

It is not difficult to verify that, if me [0, 1], then 0(m) is the unique solution of the equation 

2; = x(-T(a;),m). (8) 

At the end of the proliferating phase, a cell with a maturity m divides into two daughter cells with maturity 
g{m). We assume that g : [0,1] — > [0,1] is a continuous and strictly increasing function, continuously 
differentiable on [0, 1) and such that g{m) < m for m G [0, 1]. We also assume, for technical reason and 
without loss of generality, that 

lim g'(m) — +00. 

m — ^1 

Then we can set 

g~'^{m) — 1, for m > .g(l). 
This means that the function g^^ : [0, 1] [0, 1] is continuously differentiable and satisfies 

ig-'Yim) = 0, for m > g(l). 

Note that the maturity m of the daughter cells just after division is smaller than g{l). Then, we must 
have 

n{t, TO, 0) = 0, for TO > 5(1). (9) 
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If a daughter cell has a maturity m at birth, then the maturity of its mother at the point of cytokinesis 
was g~^{m) and, at the point of commitment, it was Q[g^^{ni)). We set 



A(m) = Q{g-^{m)), for m £ [0, 1]. 



(10) 



From a biological point of view, A gives the link between the maturity of a new born cell and the maturity 
of its mother at the point of commitment. A : [0, 1] — > [0, 1] is continuous, continuously differentiable on 
(0, 1], with A(0) = 0. Moreover, A is strictly increasing on (0,3(1)) with 8(m) < A(m) and A(m) — 0(1) 
for m G [5(1), 1]- 

Then, we can give the second boundary condition, 



n(t,TO,0) = 2(.g"i)'(m)p(i,5"i(m),r(A(m))), for t > and m G [0, 1]. 

One can note that Expression (fTT|) includes also Condition ([9|). 

To complete the description of the model, we specify initial conditions. 



p(0, TO, a) = r(TO, a), for (to, a) G [0, 1] x [0, r„ 



and 



n(0, TO, a) = n{m, a), for (to, a) G [0, 1] x [0, +C!o), 
where Tmax '■= ma-x^gjo 1] T(m) > 0. F and fi are assumed to be continuous, and the function 



II : m I 



fi{m, a)da 



(11) 

(12) 
(13) 

(14) 



is supposed to be continuous on [0, 1]. 
We put 



i{t,m) :=exp|-^ (j{x{~s,m)) + V' {x{-s,m))^ds^ , 
for t >Q and to G [0, 1], and we define the sets 

:= {(m,t) G [0,1] X [0,+(X)) ; < t < r(A(TO))|, 

and 

ne ■= |(to, t) G [0, 1] X [0, +00) ; < < < r(0(TO))|. 
Proposition 2.1. Assume that the initial conditions /i and F satisfy, for to G [0, 1], 

F(to, 0) = f3(rn,Jl{m)'j'p{m). 



(15) 



Then, the total populations of proliferating and resting cells, P{t,m) and N(t,m), satisfy, for to G [0, 1] 
and t > 0, 

d d 

—P{t,m) + —{V{m)P{t,m)) = -'y(m)P(t,m) + (3(m, N(t,m))N(t,m) 
at dm ' 



TT{m)^{t, TO )r(x(-t,TO),r(0(TO))-i), if{m,t)ene, 
7r(TO)^(r(0(TO)), to)/?(0(to), N{t ~ r(0(TO)), 0(to))) x 

iV(t-T(0(TO)),O(TO,)), ifim,t)^ne, 



(16) 
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^^N{t,m) + -^{Vim)N{t,m)) = -(^S{m) + f3{m, N{t,m))^N{t,m) 

2(.9-i)'(m)e(i, g-\m))T{x{ - t, g-'im)) , r(A(m)) - i) , tf (m, t) e Qa, (17) 

C{m)p(^A{m),N{t~T{A{m)),A{m))^N{t-T{A{m)),A{m)), if {m,t) i f^A, 

.r(e(m)) 

P(0, m) = r(m) := / r(TO, a)da, (18) 



7V(0,m) = 7Z(m), (19) 
'"^""^ " l + y(e(m))r'(e(m))' 

and 

C(m) - 2(5-i)'(m)^(T(A(m)), .g-^ (m)) . (20) 

Proof. Using (fT^ . (fT5|) and the definitions of P and N , we obtain immediatly the equations (fT5|) and (fTO]! . 
System ([I])-© can be solved by using the method of characteristics. First, we obtain the following 
representation of solutions of Equation ([1]) , 



p(t, m, a) 



^(t, m)p(0, x(— m), a — t), for < /: < a, 
^(a, m)p(t — a, x(— a, to), 0), for a < 

The initial condition (fT^ and the boundary condition give 

p(t, TO, a) = 

^(t, TO)r(x(-i,TO), a - i), for < t < a, 

(21) 

^(a, to)/3^x(— a, to), iV(i — a, x(— a, TO))^iV(i — a, %(— a, to)) , for a <t. 

Let TO e [0, 1] be given. By integrating Equation ([1]) with respect to the age, between and t(9(to)), we 
obtain 

—P(t,m) + / - — (V(m)p(t,m,a))da — —^(m)P(t,m) +p(t,m,0) - p(t,m,T(Q(m))) . 
at Jq dm ^ ' 

One can note that 

r'(e(TO))e'(TO)y(TO) - 1 = -Tx{m). 

Since 

Q rT{0{m)) Q 

— {V{m)P{t,m)) = / —{V{rn)p{t,m,a))da + T'{<d{rn))Q'{rn)V{rn)p{t,m,T{Q{m))), 

C/jn jQ 0Tn 

and 

p{t, TO, r(8(TO))) = 

' £,{t,m)T{x{-t,m),T{Q{m)) -t), if < < < r(e(TO,)), 

^(r(e(TO)), TO)/3(x(-T(e(TO)), to), 7v(i - r(e(TO)), x(-T(e(TO)), to))) x 

N{t - r(e(TO,)), x(-T(e(TO)), to)) , if T(e(TO)) < i, 
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then, using ^ and ([5]), we obtain Equation pi 
Thanks to (jlSp and by using the continuity of /I, we show that ViTHa—^+oo ^(1, m, a) = 0. 
So, by integrating Equation ^ with respect to the age, between and +oo, it follows that 

^N(t,m) + -^{V(m)N{t,m)) = - (dim) + P{m, N{t,m))]N{t,m) +n{t,m,0). 
ot dm \ J 

From the equations ([8|) and (fTO]) . we deduce that 

A(m) = x(-T(A(m)),5-i(m)). 
Hence, from Equations (fTTI) and (PT|) . we obtain 

C(to)/3(a(to), 7V(t - T(A(m)), A(m))) iV(t - r(A(m)), A(m)) , if (m, t) ^ 17a- 



n(t, TO, 0) 



Equation (fT7|) follows immediatly. 

Finally, we can remark that, if A'^ is continuous. Condition psp implies that the mappings {m,t) i— > 
F(^t,m,iV(i-T(A(TO)), A(m))^ and {m,t) G(^i, to, 7V(i - T(e(TO)), e(m)) j , with F : [0, +oo) x [0, 1] x 
R ^ M and G : [0, +oo) x [0, 1] x M ^ R given by 

f 2{g-'y{m)at,9-\m))T(x{~t,g-\m)),TiAim))~t), if(TO,t)er!A, 
F{t,m,x) = < (22) 

[ C{m)f3{A{m),x)x, if(TO,t)^17A, 

and 

( TT{m)^{t,m)r(x{ -t,m),T{e{m))-t), if (to, i) G l^e, 
Git,m,x) = l (23) 

[ 7r(TO)^(T(e(TO)), TO)/3(e(TO), x)a;, if (to, i) ^ l^e, 

are continuous. 

This completes the proof. □ 

One can remark that the solutions of Equations (fT7|) and (I19|l do not depend on the proliferating popu- 
lation. We extend N by setting 

N{t,m) ^-p{m), for t e [-r„ax,0] and TO e [0,1]. (24) 

This extension does not influence our system. However, it will be useful in the following. 

3 Local existence and global continuation 

In this section, we are interested in proving the local existence of an integrated solution of Problem (|16p - 
dH]). First, we consider an integrated formulation of Problem ([TOIl -lfT ^ . We denote by C[0, 1] the space 
of continuous functions on [0, 1], endowed with the supremum norm ||.||, defined by 

||u||= sup |w(to)|, foTveC[0,l]. (25) 

me [0,1] 

Let us consider the unbounded closed linear operator A : D{A) C C[Q, 1] C[Q, 1] defined by 
D{A) = {ue C[0, 1] ; u differentiable on (0, 1], u' e (7(0, 1], lim V{x)u'{x) = 0| 

a;— i-0 
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and 

-{5{x) + V'{x))u{x) - V{x)u'{x), if a; e (0, 1], 
-((5(0) + y'(0))u(0), ifx^O. 



Au{x) = 

Then, we have the foUowing proposition, which characterise the operator [A^ D{A)). 



Proposition 3.1. The operator A is the infinitesimal generator of the strongly continuous semigroup 
{T{t))t>o defined on C[0, 1] by 

{T{t)^lj){x) ^ K{t,x)%l;{x{-t,x)), /or e C[0,l],i > and x [0,1], 

where ^ 

K{t, x) = exp |- ^ (5{x{~s, x)) + V'{x{-s, a;)))rf.s| . 

Proof. The proof is similar to the proof of Proposition 2.4 in [TU]. □ 
We denote by C(rie) the space of continuous function on endowed with the norm ||.||f2e, defined by 

||Tl|oe:= sup |T(m,a)|, for T e C(Oe). 

(m,a)£Er2e 

Now, we can consider an integrated formulation of Problem ([TE)) - (|19p . given by the variation of constant 
formula associated to the Co-semigroup {T{t))t>o- That is the following definition. 

Definition 3.1. Let T G C(f2e) o,''^d fi be a function such that JI G C[0, 1], with Jl given by Jj^P - An 
integrated solution of Problem il6 ]) - il9\) is a continuous solution of the system 

N{t,m) = K{t,m)jl{x{-t,m)) 



Kit - s,m)/3(x(-(t - s),m),7V(s,x(-(i - s), m))) A^(s, x(-(t - s),m))ds 



(26) 



ds. 



+ K{t - s, m)F(s, xi-it - s), m), iV(s - r(A(x(-(i - s), m))), A(x(-(t - s), m)))) 
and 

P{t, m) = ^{t, m)r(x(-t, m)) 

+ I'jit ^ s,m)/3(^xi-it - s),m), Nis,xiHt - s),m))^Nis,xi-it - s),m))ds 

at - s, m)G(s, xi-it - s),m),Nis - T(e(x(-(t - s), m))), e(x(-(t - s), m))))ds. 



(27) 



for t > and m G [0, 1], where F and G are given by i2S\) and i2iJ\) and F is given by il8\) . 

The extension given by ([M)) allows the second integrals, in the expressions (pS)) and ([77]) . to be well 
defined. 

In order to obtain a result of local existence for the solutions of System ((^H|) - ([?7)) . we first focus on 
Equation ([25]) . We show, in the next theorem, that Equation ([^^ has a unique local solution, which 
depends continuously on the initial conditions. 
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Theorem 3.1. Assume that the mapping x i— > x[3{m, x) is locally Lipschitz continuous for all m G [0, 1], 
that is, for all r > 0, there exists L{r) > such that 

\xl3{m, x) — yl3{m, y)\ < L{r)\x ~ y\, if \x\ < r, \y\ < r and m e [0, 1]. 

If T E C(f2e) (ind fi is a function such that JI e C[0, 1], then, there exists Tmax > such that Equation 
i26\) has a unique continuous solution N^^'^ defined on a maximal domain [0,Tmax) x [0, 1], and either 



T 



= +00 or limsup ||iV^'^(i, Oil = +00. 



Furthermore, N^^'^{t, .) is a continuous function ofjl and T, in the sense that, ift 6 (0, Tmax), 'Pi G C'lO, 1] 
and Ti e C(rie), then there exist a continuous positive function C : [0, +oo) — > M and a constant e > 
such that, for G C[0, 1] and T2 £ C{Qe) such that A^'^2.r2 defined on [0,t] x [0, 1] and 



IIMi - /"2II < e fl"*^ llTi - r2||f2e < 

we get 

WN^^'^^is, .) - m^-'^^is, .)|| < c{t)[\\ji,-M + liri - r^llne 

Proof. We put 

r= IImII + I. 

Let T > be fixed. We consider the following set, 



forse[0,t]. 



X(7I) <^ TV G C([0,r] X [0, 1]) ; Ar(0, .) =71 on [0,1] and sup \N{t,m) --p{m)\ < 1 

[ (t,m)e[0,T]x[Oa] 

where C([0,T] x [0, 1]) is endowed with the uniform norm. X{jl) is a non-empty closed convex subset of 
C{[0,T] X [0,1]). 

We define the operator H : C([0, T] x [0, 1]) ^ C([0, T] x [0, 1]) by 
H{N){t, m) - K{t, m)-p{x{~t, m)) 

K{t - s, m)(3(x{-'{t " s), m), 7V(s, - s), m))) iV(s, - s), m))ds 

+ K{t - s, m)F(s, x{-{t - s), m), N{s - r(A(x(-(i - s), m))), A(x(-(i - s), m)))) ds. 

H is continuous in C([0,r] x [0,1]). Our objective is to show that iJ is a contraction from X{'fl) into 
itself. 

Let N e X{ji). It is clear that H{N){0, .) = /I. On the other hand, we have, for {t, m) e [0, T] x [0, 1], 
\H{N){t,m) —Ji{m)\ < | TO)7t(x(^i5 "i)) ^ Pifn)] 



K 



{t - s, m)/3(x(-(t - s), m), 7V(s, - s), m))) Ar(s, x(-(t - s), m))ds 

- s, m)F(s, - s), m), iV(s - T(A(x(-(t - s), m))), A(x(-(t - s), m)))) 



ds 
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Since K is continuous on [0, T] x [0, 1], then there exists K >Q such that 

\K{t,m)\< K, for (t, m) G [0, T] x [0, 1] . 

Since N e X{ji), then 

\N{t,m)\ < 1 + \\J1\\ = r. 

This imphes that there exists M := max{M, ||C||rL(r)} > such that, for {t,m) e [0, T] x [0,1] and 

se [o,t]. 



F 



(s, x{-(t - s), m), N{s - r(A(x(-(t - s), m))), A(x(-(t - s), m)))) 



< M, 



where 



M := sup 

(m,t)fEnA 



and C is given by ([^0]) . Hence, we obtain that 

|i7(iV)(i,m) - 7I(to)| < \K{t,m)jl{x{~t,m)) -Jl{m)\ + K{rL{r) + M)t. 

Let us recall that K{0,m) = 1, x(0, m) = m and {t,m) t— K{t,m)'Jl(^x{~tjiT^)) is continuous. Then, we 
can choose T > such that 



Consequently, 



sup < \K{t,m)jL{x{'-t,m)) - 'p{m)\ + KlrLlr) + M)t\ < 1. 
(t,m)e[o,r]x[o,i] 



\H{N){t,m) ~Jl{m)\ < 1, for {t,m) G [0, T] x [0,1], 



(28) 



and H{X{-p)) cX(p). 

Now, we show that iJ is a contraction on X{jl). 
Let Ni e XijL) and A^a S X(7l). Then, 

|i7(7Vi)(i,m)-i7(A^2)(t,m)| 



< 



iir(i — s, m) 



/3(x(-(t - s), m), iVi(s, x(-(t - s), m))) Ari(s, x(-(t - s), m)) 
-P{x{~{t - s), m), iV2(s, x{-{t - s), m)))iV2(s, x(-(t - s), m)) 



— s, m) 



(s, x(-(t - s), m), 7Vi(s - T(A(x(-(i - s), m))), A(x(-(t - s), m)))) 



x{-{t - s), m), 7V2(s - r(A(x(-(t - s), m))), A(x(-(t - s), m)))) 



< is:(l + ||C||)i(r-)T sup \Ni{t,m)-N2{t,m)\. 

(t,m)e[0,T]x[04] 

Since r > 1 and ||C||rL(r) < M, then Condition implies that 

^(1 + ||C||)i(r)T < ^(1 + ||C||)rL(r)T < K{rL{r) + M)T < 1. 



10 



M. Adimy and F. Crauste 



A nonlinear model of cellular proliferation 



Hence, H is a contraction from X{ii) into itself. Therefore, there exists a unique N G X{ii) such that N 
satisfies Equation ^ for (t, m) e [0,T] x [0,1]. 

Let N be the maximal solution of Equation (|26p . defined on its maximal domain [0,T„iq2;) x [0, 1]. 
Assume that 

Tmax < +00 and limsup ||A^(i, .)|| < +oo. 
Then, there exists r > such that 

||iV(i,.)|| <r, fort e [0,T„,,). 
Let {t, m) e [0, T^ax) x [0, 1] and h>{) such that t + e [0, Ty^ax)- Then, 
\N{t + h,'m) - N{t,m)\ 

< K{t + h, m)Ti{x{-{t + h),m)) - K{t, m)Ti{xht, m)) 
ft+h 

K{t + h - s,m)f3lx{-{t + h- s), to), 7V(s, x(-(t + h- s), to)) j iV(s, x(-(t + h- s),m))ds 
K{t - s, m)/3(x(-(t - s),m), N{s, x{~{t - s), TO)))iV(s, x{~{t - s), m))c 

t+h 

K{t + h — s,to)x 



^(s, Xi~{t + h-s), m), N{s - r(A(x(-(< + - s), m))), A(x(-(t + - s), m))))ds 
if(t - s,m)F(s,x{-{t - 5),to),7V(5 - r(A(x(-(t - s) , m))) , A{x{~ {t - s),m))))ds 



< 



K{t + h, m)jL{x{-{t + h), to)) - K{t, m)jl{xi-t, m)) 
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Since 



h := K{t + h- s, m)/3(x(-(i + h-s), m), N{s, x{-{t + h-s), m)))iV(s, x{-{t + h-s), m))ds 

-J^K{t- s, m)/3(x(-(t - s), m), N{s, x(-(t - s), m))) Ar(s, x{-{t - s), m))ds, 

= J ^K{t-s, m)/?(x(-(t - s), m), N{s + h, x(-(t - s), m))^N{s + h, x(-(t - s), m))ds 

- j^K{t- s, m)(i(x{-{t - s), m), N{s, x(-(i - s),m)))N{s, x(-(i - s), m))ds, 

= J ^K{t-s, m)/3(x(-(i - s), m), N{s + h, x{-{t - s), m))) Ar(s + h, x(-(i - s), m))ds 



+ / K{t- s,m) 
Jo 



p{x{-{t - s), m), N{s + h, x{-{t - s), m)))iV(s + h, x{-{t - s), m)) 
-P{x{-{t - s), m), N{s, x(-(i - s),m))')N{s, x(-(t - s), m)) 



then 



ds. 



\h\<KrL{r)h + KL{r) / N {s + h,x{-{t - s),m)) - N {s,x{-{t - s),m)) 

Jo 



ds. 



By using the same reasonning, we also obtain that 



1/2 1 < 



KMh + giih) + K\\C\\L{r) f iV(s + /i - r(A(x(-(t - s),m))), A(x(-(i - s),m))) 

Jo 



-N{s - r(A(x(-(t - s),m))), A(x(-(t - s),m))) 



ds. 



with gi{h) independent of (t, m) and such that Hm^i^o 8i{h) — 0. 
One can remark that, if ti > 0, ^2 > and m € [0, 1], then 

K{ti + 12, m) = K{h,m)K{t2, xHum)) = K{t2, m)K{ti,x{-t2, m)). 

This yields to 

\N{t + h,m) - N{t,m)\ 

< K(Q2{h) + rL{r)h + Mh) + Q^{h) + kL{r) \\N{s + h,.)- N{s, .)\\ds 
+K\\C\\L{r) ||iv(s + h- r(A(x(-(t - s), m))), .) - 7v(s - r(A(x(-(t - s), m))), .) 



ds, 



where 



Q2{h) = sup 

(t,m)e[0,T„„i,)x[0,l] 



K{h, xi-t, m))-p{x{-{t + h),m)) - -p{x{-t, m)) 



h^O 
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We set _ 

g{h) = K(^Q2{h) + rL{r)h + Mh^ + gi{h). 

Hence, for {t, m) e [0, Tmax) x [0, 1], > and 9 e [-T,nax, 0] such that t + he[0, T^ax) and t + 9>0, 

\N{t + h + e,m)- N{t + e,m)\< g{h) + KL{r) / \\N{s + h, .) - N{s, .)\\ds 

Jo 

j-t+e 

+K\\aL{r) / \\N{s + h- r(A(x(-(< - s),m))), .) - N{s - r(A(x(-(< - s),m))), .) 
Jo 

On the other hand, 

N{t,m) —Jl{m), for t G [— Tv„Qa;,0] and m e [0,1]. 

This iniphes that 

sup \\N{t+h+9,.)-N{t+e,.)\\ < g{h)+K{\\C\\ + l)L{r) f sup \\N{s+h+0, .)~N{s+0, .)\\ds. 

By using the Gronwall's inequahty, it foUows that 

sup \\N{t + h + 9,.)-N{t + d, .) II < p(/x)e-'f (IICII+i)^^^™.^ _ 
ee[-r,„„,,o] 

Hence, 

||iV(i + /i, .) - N{t, Oil < e(/i)e^(llCll+i)^('-)^"- , 

with 

hm g(h) = 0. 

Using the same reasonning, we can show a same result for ft, < 0. 
It foUows immediatly that 

hm N{t, .) exists. 

This imphes that N can be extended continuously to Tmax, which contradicts the maximality of [0, Tmax)- 
Hence, 

Tmax — ~^oo or lim sup 1 1 A^(t , .) || = +CXD. 

Now, we denote by N^^^'^^ the solution of Equation (|26|) for the initial data /Ii and Fi, defined on its 
maximal domain \}i,Tmax) x [0, 1]. Let t G {Q,Tmax) be fixed. We put 

r(i) = sup ||iV'^i'^i(s,.)|| and R{t)^l + r{t). 

se[o,t] 

We set _ 

C:= sup |2(5-i)'(m)e(t,.g-i(m))|, 

and 

C[t) = Xmax{l,C}e^(«('»^(i+ll';il)*. 

Let < e < 1 be such that 2C{t)e e (0, 1) and let Ji^ e C[0, 1] and Fa G C(f2e) such that Af'^a-rs is 
defined on [0, t] x [0, 1] and 

IIMi-M2II<£ and HFi - FaHne < £• 
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Then, there exists s > such that ||iV^2,r2(^^ )|| < ^(^) foj. a e [0, s]. Let 
to = sup|s>0; ||7V'^2'^^(ct, Oil < for cr e [0, s] 

If we suppose that to < t, then N^^'^^ is defined on [0, to] x [0, 1] and satisfies, for s G [0, to] and m G [0, 1], 

\mi.^^{s,m)-m^^^^{s,m)\ < K(Wjl,-Jl2\\+C\\T,-T2\\ne) 

+ Kil + \\C\\)LiRit)) r\\N^^-^^{a,.)~N^^-^Ha,.)\\da. 
Jo 

Therefore, using the GronwaU's Inequahty, we obtain 

\m^^^'is,rn) - m-'^-{s,rn)\ <C{t)(^\\jl,~Jl2\\ + \\Ti-T2\\no), forsG [0,io]. (29) 
This imphes, in particular, that 

\\N^''^Hs,.)\\ <2Cit)e+\\m^-^'is,.)\\ < l + r(t) for s G [0,io]- 

This contradicts the definition of to- Hence, to > t. 

That means that ((^ is satisfied for each t G [0,Tmax)- Then, we deduce the continuous dependence of 
the solution with the initial data and the proof is complete. □ 



One can remark that Condition (fT5|) is not needed to prove Theorem 13. II 

By using Theorem 13. 11 we can deduce the following result, which deals with the existence of solutions 
of Problem 

Corollary 3.1. Under the assumptions of Theorem \3.1\ System i26]) - (W^ has a unique continuous max- 
imal solution {N,P), defined on [0,Tmax) x [0,1]- 

Proof. Under the assumptions of Theorem l3.1[ Equation (j26p has a unique continuous maximal solution N 
defined on [0, T„iax) x [0, 1]. Then, we easily obtain the existence and uniqueness of a solution of Equation 
(EZl) on [0,T^a,) X [0,1]. □ 

We can use the results of local existence, given by Theorem 13.11 and CoroUarv 13.11 to investigate the 
global existence of the solutions of System ([26 ]) - ([27]) . This is done in the next theorem. 



Theorem 3.2. Under the assumptions of Theorem \3.1\ and the assumption that the mapping x t— > /3(m, x) 
is uniformly bounded, the unique solution of System Il26\ )- (2l\ l is global, that means, it is defined for all 
t > 0. 

Proof. We assume that there exists /? > such that 

]P{m,x)\ < /3, for aU m G [0, 1] and x G M. 
Then, for t G [0, Tmax) and m G [0, 1], we get 

\N{t,m)\ < KWn+Kp f' \\N{s,.)\\ds 







+ KMt + K(3\\C\\ / \\N{s - T{Aixi-{t - s),m))), .)\\ds. 
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Let 6 e [-Tmaa;, 0] be such that t + 9 >Q. Then, 

\N{t + e,m)\<K{\\TA\+Mt)^kp{\\(:\\ + l) [ _ sup \\N{s+'e,.)\\ds. 

On the other hand, if —Tmax < t + 9 < 0, then 

\N{t + 9,m)\<Km\, 

because K > 1. Hence, 

sup \\Nit + 9,.)\\<Km\+Mt)+KM\\C\\ + l) f sup \\Nis + 9, .)\\ds. 
0e[-T,„ax,o] Jo ee[-T,„a^,o] 

By using the Gronwall's inequahty, we deduce that 

\\N{t,.)\\< sup ||iV(i + 0,.)|| <i^(||M||+Mi)e^^(ll«ll+i)* :=/(t). 

Since / is continuous on [O^Tmax], then 

hmsup \\N{t, .)|| < +00. 
We deduce that T^ax — +oo and the solution of (|26p is globaL 

FinaUy, we easily obtain that the unique maximal solution of Equation (f27|l is also global. □ 

Corollary 3.2. Assume that the mapping x i— > x(3{m, x) is Lipschitz continuous for all m £ [0, 1]. Then, 
for r e C(r2e) and fj, such that JI G C[0, 1], the unique solution of System i26]) - (W^ is global. 

We have studied, throughout Theorems 13.11 and 13 . 2[ the local and global existence of the solutions of 
((25)) and of Problem (P^ - ([?7)) . Before we investigate the positivity of these solutions in the next section, 
we can ask for regularity results for these solutions. This is presented in the following remark. 

Remark 1. Under classical assumptions on the function (3 and the initial data T and /i, we can obtain 
regularity results for the solution of System i26] )- (W^ . This may be done by using the same idea as in the 
work of Travis and Webb J23^ . 

4 Positivity of solutions 

Since we study a biological population, it is necessary to obtain the positivity of the solutions of System 
([26| - ([27|) to ensure that the model is well-posed. First, we focus our study on the solutions of Equation 
([26| . We use a method given by Webb [24] in 1985 and developped by Kato [14 to obtain the positivity 
of the solutions of Equation ([25]) . 

Let r e C(rie) and /i be a function such that JI G C[0, 1], JI given by (HH). Let T > be fixed. 
We consider the family of operators H'' : C([0,T) x [0,1]) C([0,T) x [0,1]) defined, for a e C[0, 1], 
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{t,m) G [0,T) X [0, 1] and N e C([0,T) x [0, 1]), by 

H''{N){t,m) = e-*'»W-*'™))ii'(i,m)7Z(x(-i,m)) 



e-(*-«)«W-*'™»/f(t-s,m) a(x{-{t-s),m))N(^s,x{-{t-s),m)) 
-^(xi -{t- s),m),7v(s,x( - (i - s),m))^7v(s,x( - {t - s),m)) 

X( - - s), m) , 7V(s - T(A(x(-(f - s), m))), A(x(-(t - s), m)))) 



(30) 



Let a e C[0,1], (t,m) e (0,T) x [0,1] and iV e C([0,T) x [0,1]) be fixed. We consider the mapping 
w" : [0,t] ^ K defined, for s G [0,t], by 

w-{s) = H-{N)(s,xi-{t-s),m)). 

Tlien, we can prove the following lemmas. 

Lemma 4.1. The function w" is differentiable on [0,t] and satisfies 

( d 



ds 



w"(.s) = -!^"(s)u;"(s) +/"(s), forse[0,t], 
'^'"(O) = -p{x{-t,m)), 



with 



and 



u'^is) = a{x{-t,m)) + S(^x{ - {t - s),m)) + (x( ~ {t - s),m)). 
Proof. First, one can remark that, if ti and t2 arc positive, then 

x( - (ti+t2),m) = x{-ti,x{-t2,m)). 
It follows that, for s G [0,t] and h^O such that s + h G [0,t],we have 

1 



with 



h 



(31) 



r{s) = a{x{ - it - s),m))N(^s,x{ - {t - s),m)') 

- n(^s, x{-{t- s), m))/?(^x( -{t- s), m) , 7v(s, x( - (i - s), m)) 

+ F^s,x( - (t - s),m), Ar(^s - t(a(x( - (i - s),m))), a(x( - (t - s),m))^ j . 



+ /i) - w^"(s)) = Wi,h + W2,h, 



ha(x{-t,m))j^(^^^^^ - (t - s - h),m)) - 1 e-««W-*'™))i^(s,x( - - s),m))7l(x(-<,m)). 
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and 



W- 



2,h 



= - J e-("+'*-^)"(^(-*''"»if(s + /i-CT,x(-(t-s-/i),m)jx 
a(^xi-{t-'j),m)^N(a,x{-{t-a),m)^ 
-N(a,x{ - (i - (t), m))/3(^x( - {t - (t), m) ,N(a, x{-{t- a), m)) 
+F^a,x{-{t-a),m),N(^a-T(^A{x{-{t-(7), 



m))),A{x{-{t-a),m 



m X 



a(x( - (i - <^), m)^N(a, %(-(*- CT),m)) 
-iv(a,x(-(i-^^),m))/3(^x(-(i-^^),H'^(^'X(-(i-«^),m))^ 

+F (^a,x{ -{t- a), m),N(^- t(a(x( - (i - a), m))) , a(x( - (t - a), m 

Hence, we easily obtain that 

lim W^i,;, = -zv«(s)e--(>^(-*'™»if(s,x( - (i - 5),m))M(x(-t,m)). 



da 



da. 
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By the same way, we have 
1 



Jo 



h 



-(s-<7)a(x(-t, 



»i^(s-a,x(-(i-s),m)) a(x( - (t - a), m))7v(a, x( - - a), m)) 



-N(a,x{ - {t - a),m)')0(x{ - {t - a),m),N(a,x{ - {t - a),m))i^ 

+F ^a, X ( - (t - a) , m) , TV (^a - T ( A (x ( - (i - , m) ) ) , A (x ( - (t - ^) , m) ) ^ j 



da 



+ ^e-'*°W-*''"»i^(/i,x( - (i - s - h),m)^ ^ e-^s-a)a(x(-t,m)) _cr,x{-{t- s), 

a{x{ -{t- (t), m))iv(a,x( - {t - a),m)) 
-n(c7, x{-{t- (t), m))/3(^x( -{t- a), m) , iv(a, x( - {t - a), m))^ 
+F ( a, X ( - (i - a) , m) , A/- (^a - r ( A (x ( - (i - a) , m) ) ) , A (x ( - (t - a) , m) ) ^ 



m X 



Therefore 



lim W2,/i 

?i->0 



da. 



-z/»(s) y e-(^-'^)"(''(-*'"'»Js:(s-(7,x(-(i-s),m))x 

a (x ( - (* - , ) AT (<T, X ( - (i - ct) , m) ) 
-N(a,x{ - {t - a),m)^p(^x{ - {t - (r),m),N(a,x{ - {t - a),m)^^ 
+F l^a, x( - (i - a), m), iV (^ti - t(a(x( - (i - a), m))) , a(x( - (t - d), m))^ 
+a(x( - (i - s), m)) 7v(s, x( - (< - s), m)) 

-iv(s, x( - (i - s), m))/3 (^x( - - s), m) , A^(s, x( - - s), m))^ 

+fIs,x{ -{t- s),m),7v(^s - t(a(x( - {t - s),m))), a(x( - - s),m))^ j 



da 
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Hence, lim^j^o ^(^'^(^ + h) — w°-{s)) exists and is equal to — i^"(s)w"(s) + /"(s). 
We obtain Equation (|3ip and the proof of the lemma is complete. 



□ 



Lemma 4.2. Let a,b E C[0, 1]. Then, 

H\N){t,m) ^ H''{N){t,m)+ f e-(*-")''('^(-*'"'»i^(i - s, m) x 

(&(x(-(t-.s),m)) -a(x(-(t-s),m))) iv(s, x( - (t - s), m)) ~w^{s) 



ds. 



Proof. From Equation ([31]) . it follows that 
d 



ds 



w 



+ {b{xi-{t - s),m)) ^ a{xi-{t ^ s),m))^ n(^s,x{~ (t - s),m)^ 

Since w''{0) = w"-{0) = Jl{xi-t, m)), then 
{w''it) -w''{t))e^o''''('^')da ^ 

Jo '^'(-)d- (b{x{-{t - s),m)) - a{x{-{t - s), m))) x( - (* - s), m)) ~ 



ds. 



Therefore, 

w''it) -w^it) = 

This ends the proof. 



e-it-'^Wx{~t^"^))K{t - s,m)x 
{b{x{~{t - s), m)) - a{x{-{t - s), m))) iv(s, x( - (i - s), m)) - w''{s) 



ds. 



□ 



Lemma 4.3. Assume that there exist a G C[0, 1] awd N e C([0,T) x [0, 1]) such that H°-{N) = N. Then, 
for all b e C[0, 1], H''{N) ^ N. 



Proof. This is a simple consequence of Lemma 



□ 



Now, we show, under a cellular regulation hypothesis, that the solutions of (P5|l -(|27 p are positive. We 
use the three previous lemmas to show the following theorem. 

Theorem 4.1. Assume that the mapping x ^ f3{m,x) satisfies the cellular regulation hypothesis, 

{/3{m, x) - (3{m, 0))x < 0, for m e [0, 1] and x G R. (32) 

Let r e C(r2e) and ^ he a continuous function such that JI G C[0, 1]. // /I > and F > 0, then the 
solutions of Equation i26]) are non-negative on their domains. 

Proof Let TV be a solution of Equation ^ defined on [0,T) x [0, 1], where T > 0. 
We consider the operator defined, for a G C[0, 1], by ([50]) . It is obvious that 

H"{N) = N, 



where denotes the null mapping of C[0, 1]. 
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Then, from Lemma for all a e C[0, 1], H°-{N) = N. Hence, if there exists a e C[0, 1] such that 
H''{N)it,m) > Ofor {t,m) e [0,r)x [0, 1], then we shaU obtain that N{t,m) > 0, for (t, m) G [0,T)x [0, 1]. 
We set 

a(m) := /3(m, 0), for aU m £ [0, 1]. 
Then, a G C[0, 1] and Condition ((32|) implies that 

a{m)N{t,m) - N{t,m)/3{m, N{t,m)) > 0, for (t,m) G [0, T) x [0,1]. 

Hence, 



a[x{ -{t- s), m))iv(s, x( - - s), m)) 



-/3(^X( - - s), H' ^(^' x{-{t' s), m))^iv(s, x( - (< - s), m)) 
is non-negative. We have to show by steps that, for {t, m) G [0, T) x [0, 1], 

g-(t-.)a(x(-t,m))^(^^ _ ^^p(^g^ x{ -{t- s), m),N{s - T{A{x{-{t - s),m))),A{x{~{t - s), m))))ds 



is non-negative. 

First, it is clear that H°'{N) is non- negative on J^a- 
We set 

TA := min rfAfm)). 

me[0,l] 

Since t is positive on [0, 1], then ta > 0. 

Let m G [0, 1]. Suppose that t G [r(A(m)), r(A(m)) + ta]. Then, 

< s - t{A{x{ - (< - s),m))) < t - TA < r(A(m)), for s G [T(A(m)),t]. 
Since _F > for (m,t) G JIa, this yields that 



t 

e' 



(t-.)a(x(-t,m))^^^_^^ m)F(s, x(-(t-s), m), iV(s-T(A(x(-(t-s), m))), A(x(-(t-s), m))))ds > 0. 



Hence, H''{N){t, m) > for m G [0, 1] and t G [0, T{A{m)) + ta]. This implies that 
N{t, m) > 0, for m G [0, 1] and t G [0, r(A(m)) + ta]. 

By steps, we show that 

N{t, m) > 0, for m G [0, 1] and t G [0, t(A(to)) + hta] with n G N. 
There exists n G N such that 

r(A(m)) + nTA < T < t(A(to)) + {n + 1)ta. 
Consequently, > on its domain [0, T) x [0, 1] and the proof is complete. □ 
We deduce immediatly the following corollary. 



Corollary 4.1. Under the assumptions of Theorem \4-l\ the solutions of Problem i26]) -M^'AI are non- 
negative. 
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5 Local and global stability 

We study, in this section, the local and global stability of the trivial solution of Problem (P5)l - (|27p . We 
first focus on Equation (|26p . 
We set 

i5 := inf \5{m) + V'{m)], k := sup \{g~^)'{m)\ and 7:= inf [^{ra) + V'{m)]. 
me[0,l] V / me[0,l] me[o,i] V / 

We assume that 7 > and that the function x ^ xl3{m, x) satisfies a Lipshitz condition in a neighborhood 
of zero, that is, there exist e > and L > such that 

\x(3{m,x) — y(3{m,y)\ < L\x — y\, if < e, |y| < e and m G [0, 1]. (33) 

Throughout this section, N^^^^ denotes a global solution of Equation ([26|) associated with the initial data 
^ and r. 

In order to obtain the local stability of the trivial solution, we first prove a result of invariance for the 
solutions. 

Proposition 5.1. Assume that Condition hS^] is satisfied and that 

L{1 + 2k) < S. (34) 

Let r G C{fls) arid ^ be a continuous function such that JI G C[0, 1]. // ||7l|| < £ and ||r||sig, < eL, where 
e and L are given by i33\) . then 

|A^^'^(t,m)| < £, fort>0 and m G [0, 1]. 

Proof. We first notice that, under the assumption (|34|) . (5 > 0. We consider the sequence defined by 

No{t, m) = K{t, m)jl{x{-t, m)) , 

and 

Nnit, m) ^ No{t, m) - /(iV„_i)(t, m) + J(7V„_i)(i, m), 
for t > 0, TO G [0, 1] and n eW, with 

I{N)it, m)^ Kit - s, m)f3[x{-it - s),m),N{s, x(-(i - s), to))) A^(s, x(-(i - s), m))ds, 

and 

JiN){t, m) = ^ Kit - s, m)F(^s, xi-{t - s),m),N{s - r(A(x(-(i - s), to))), A(x(-(t - s), TO))))ds. 
Then, 

\Noit,m)\ < e"**||7i|| < e-^^e < e. 

By induction, we show that 

\Nnit,m)\<e, for t > 0, TO G [0,1] and n G N*. (35) 

We assume that |iV„(<, to)| < e. 

Let it, to) G [0, +00) X [0, 1] be fixed. 
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One has to notice that, if < s < t(A(x(— (i — s),r7i))), then 

\F(s,x{-{t - s),m),7V„(s - r(A(x(-(t - s),m))), A{x{~{t - s),m)))) \ 
<2ne-nnne, 

and, if s > t(A(x(— (i — s),r7i))), then 

\f(s, xi-{t - s), m), N^s - r(A(x(-(< - s), m))), A(x(-(t - s), m)))) | 

< IICII /3(A(x(-(t - s), m)), iV„ (,s - r(A(x(-(t - s), m))), A(x(-(i - s), m)))) : 



iV„(s - r(A(x(-(t - s),m))), A(x(-(t - s),m))) 
< KWeL, 

where ||C|| is given by ([^5]) . Since \\(\\ < 2k, we obtain that 



|J(iV„)(i,m)| <2KeL / e-^'^'-'^ds, 



and 



|iV„H.i(t,m)| < \No{t,m)\ + \I{N„)(t,m)\ + \J{Nr,){t,m)\, 
< ee-^* + ei(l + 2K) / e~^'^*~^Us, 



< ee-** + er^(l + 2K)(l-e~**) 
(5 



By assumption. 
Consequently, 



L 



(1 + 2k) < 1. 



\Nn+iit,m)\ <£(e-'^*(l-£(l + 2K)) +£(1 + 2k)) < £. 



We conclude that is true for n € N* . 

By remarking that the sequence {Nn)nefi* converges to N^'^ , then we obtain that 



|7V^'^(t,TO)| < e, for t > and m G [0, 1] 



which ends the proof. 



□ 



Before we prove the local stability of Equation , we rewrite the integrated solution of by using 
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a variation of constant formula, for t > Tmax and m G [0, 1]. Wc obtain that 

N{t,m) ^ K{t - Tmax,m)N(Tmax,x{ " " Tmax) , m)^ 

- ^ K{t- s, m)/3(x(-(t - s), m), N{s, x{-{t - s), m))^N{s, x{~{t ~ rn))ds 

+ f K{t- s, m)C(x(-(i - s), m))N{s - r(A(x(-(t - s), m))), A(x(-(t - s), m))) x 

p(A{x{-it - s),m)),N{s - r(A(x(-(t - s) , m))) , A{x{~ {t - s), m))))ds. 

Hence, we can show the foUowing theorem, which deals with the local stability of the trivial solution of 
Equation (|26l) . 

Theorem 5.1. Under the assumptions of Proposition [KT\ the trivial solution of Equation 126]) is locally 
exponentially stable, that is, there exist t > 0, e > 0, c > and d > such that, if \\jl\\ < e and ||r||ne < 
then 

WNf'-'^ {t, .)\\ < ce-'^'^'-~'\ fort>t. (36) 
Proof. Let e > and L > be given by ([33|) . We assume that 

||7I|| < e and ||r||oe < eL. 

Then, Proposition 15. II implies that 

\N'>^^^{t,m)\ < e, for i > and TO e [0, 1]. 
We define the sequence (iV„)„gN by 

NQ{t,m) = K{t ~ Tmax,m)N'''^ (Trnax,x{ ^ ^ Tmax) , m)^ , 

and 

Nn{t, to) = No{t, m) - (Nn-i){t, to) + J^„_ (Nn-i){t, m), 
for t > Tmax, m £ [0, 1] and n eW, with 

iN)it, m)^ j Kit - s, m)/?(x(-(i - s), m),N{s, x(-(t - s), to))) 7V(s, x(-(i - s), m))ds, 

and 

Jr^.AN)it,m) - / K{t-s,m)ax{-{t-s),m))N{s~T{A{xi-{t-s),m))),A{xi~{t-s),m)))x 

J Tmax 

P[A{x{-{t - s),m)),N{s - r(A(x(-(^ - s), to))), A{x{-{t - s),mj))ys, 

and 

Nn{t, to) N^-'^{t, to), for [t, to) e [0, Tmax] X [0, 1] and n e N. 

First, it is easy to check, by using the same reasonning as in the proof of Proposition IS . 1 [ that the sequence 
{Nn)neti satisfies 

\Nn{t, to)| < e, for t > 0, TO e [0, 1] and n e N*, 
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since 

\No{t, m)\<e, for t > and m e [0, 1]. 

Secondly, the assumption ^ 

L(1 + 2k) < 6, 

implies that _ 

Lil + \\C\\)<5. 

Therefore, there exists p £ (0, d) such that 

S-p S 

L < TTTTT- < 



l + ||^||ePWx i + ||^||- 
Let t > Tmax and m G [0, 1]. Then, wc get 

\Noit,m)\ < e-^(*-^— )£ < e-''(*-^-'"')e < s. 

Since 

\N,{t,m)-No{t,m)\ < |/r_.(^o)(i,m)| + |J._(iVo)(i,m)|, 

then, the estimates 

\Ir^a. (No){t, m)\<eL t e-'^'-^^e-^^^-^-^Us, 

and 

t 



yield to 



\Ni{t,m) - Noit,m)\ < eL{l + HCHe"^— )e-''*e''^— / e^^'P'^'ds. 
By remarking that 

J Tmax 5 — p 

we obtain 

\Ni{t,m)-No{t,m)\ < eL^-^W^—^e-''^*-^'-^-\ 

S-p 

It is easy to see by induction, that 

|7V„(t,m) - iV„_i(t,m)| < ej^Lil^lif^^^ e-p(t-r„ax)^ for ^eW 

Hence, we get 

n 

\Nn{t,m)\ < Y^\Ni{t,m)-Ni-i{t,m)\ + \No{t,m)\, 

1 +11(116"^- 



i=0 



S-p 
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Since 



then, we finally obtain that 



|7V^-r(^ )| < ^ ^-P{t-r^..) 

Hence, we have obtained the inequality p6p with t = Tmax, e = niax{e, eL}, d = p and 

^-p-Lll + llClle"-"-)' 

This completes the proof. □ 

By using Theorem 15.11 we can prove the local exponential stability of the trivial solution of Problem 
((?f|) . This result is presented in the following corollary. 



Corollary 5.1. Assume that 

L(l + 2k) < min{7, (5}. 
Then, the trivial solution of Problem i26] )- (W^ is locally exponentially stable. 

Proof. First, since L(l + 2k) < S, then Theorem 15.11 implies that the trivial solution of Equation (|26p is 
locally exponentially stable. Moreover, since we suppose that L(l + 2k) < 7, we can show, by using the 
same arguments as in the proof of Theorem 15.11 and a variation of constant formula for P, that the trivial 
solution of Equation (I7f|) is locally exponentially stable, as soon as the trivial solution of Equation 
has this property. □ 



Finally, by using Theorem 15.11 we can obtain the global stability of the trivial solution of Problem 
26 | - (|27|) . This is done in the next proposition. 



Proposition 5.2. Assume that the mapping x ^ xP{m,x) is Lipschitz continuous for all m e [0, 1], with 
a Lipschitz constant L, and that 

L{1 + 2k) < mm{j,5}. 

Let r £ C(Qq) and fi be a continuous function such that JI £ C[0, 1]. // ||r||s7g, < then, the trivial 

solution of Problem i26\} - (27\ l is globally exponentially stable. 



Proof. We first show, as in the proof of Proposition [531 that 

|iV'^'^(i,m)| < for t > and m G [0,1]. 

Then, by using the same reasonning as in the proofs of Theorem 15. II and Corollarv l5.ll we conclude. □ 

We can give some explanations about the condition i(l + 2k) < mhi{j,d}. In fact, we can notice 
that this inequality is satisfied if 7 and S are large or if L is small enough. This corresponds, biologically, 
to the case where the mortality rates (7 and d) are important or to the case when only a few cells are 
introduced in the proliferating phase, and then the cells supply is not sufhcicnt {L is a bound of the 
number of introduced cells). 
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6 Discussion 

It is usually believed that the function /3 is a Hill function (see Mackey [THj) given by 

Mm) ^ , , — , for X > 0, m e [0, 1], 
/3o(to), for X < 0, m e [0, 1], 

with and Pq two continuous and positive functions, and n > 1. In this case, the mapping x xf3{m,x) 
is always Lipshitz continuous, with a Lipschitz constant equals to 

sup /3o(to). 

m6[0,l] 

Hence, the results of Theorem 13. II and Corollarv l3.2l hold. 

By the same way. Condition (l32|) is easily satisfied, because /3(to, x) — (3{m, 0) if x < 0, and, if x > 0, 
the function x i— > (3{m,x) is decreasing for all m. Hence, the positivity of the solutions is naturally 
obtained. 

As it has already been noticed by Dyson et al. [TU], Mackey and Rudnicki [2D] and Adimy and Pujo- 
Menjouet [3], we can expect to show the influence of the immature cells population (that means, the 
population with a small maturity) over the entire population. This has also been obtained by Adimy and 
Crauste [1] for a model with a proliferating phase duration distributed according to a density. 

In particular, by using the result of local stability obtained in Section O we will certainly be able to 
prove that the stability (or the instability) of the immature cells population leads to the global stability 
(or the instability) of the entire population. This has been displayed for the first time by Mackey and 
Rudnicki [501 i'^ 1999. Since the stem cells population is known to be at the root of the blood production 
system, then, this behaviour is naturally expected in our model, and it is the purpose of a next work. 
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